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ABSTRACT 

We assess a model of late cosmic reionization in which the ionizing background radiation arises 
entirely from high redshift quasars and other active galactic nuclei (AGNs). The low optical depth 
to Thomson scattering reported by the Planck Collaboration pushes the redshift of instantaneous 
reionization down to z = 8.8+14 and greatly reduces the need for significant Lyman-continuum emis¬ 
sion at very early times. We show that, if recent claims of a numerous population of faint AGNs at 
z = 4 — 6 are upheld, and the high inferred AGN comoving emissivity at these epochs persists to 
higher, z > 10, redshifts, then active galaxies may drive the reionization of hydrogen and helium with 
little contribution from normal star-forming galaxies. We discuss an AGN-dominated scenario that 
satisfies a number of observational constraints: the H 1 photoionization rate is relatively flat over the 
range 2 < z < 5, hydrogen gets fully reionized by z ~ 5.7, and the integrated Thomson scattering 
optical depth is r ~ 0.056, in agreement with measurements based on the Lya opacity of the inter- 
galactic medium (IGM) and cosmic microwave background (CMB) polarization. It is a prediction of 
the model that helium gets doubly reionized before redshift 4, the heat input from helium reionization 
dominates the thermal balance of the IGM after hydrogen reionization, and z > 5 AGNs provide a 
significant fraction of the unresolved X-ray background at 2 keV. Singly- and doubly-ionized helium 
contribute about 13% to r, and the He hi volume fraction is already 50% when hydrogen becomes 
fully reionized. 

Subject headings: cosmology: theory — dark ages, reionization, first stars — diffuse radiation 
intergalactic medium — galaxies: active 


1. INTRODUCTION 

The reionization of the all-pervading IGM marks a 
turning point in the history of structure for nation in 
the universe. The details of this process refle ;t the na¬ 
ture of the first astrophysical sources of radiition and 
heating as well as the early thermodynamics of cosmic 
baryons, and continue to be the subject of considerable 
obs ervational and theoreti cal efforts (for a rece at review, 
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the integrated UV emission from active nuclei and star¬ 
forming galaxies, but there is still no consens us on the 
relat ive contributions of these sources as a fi motion of 


epoc h The high ionization threshold (4 ryd) 411 d small 
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photoionization cross section of He 11 , the ra; 
bination rate of He hi, and the fact that most 
lack 4 ryd emission, all delay the double ioni: 
helium. This is expected to be completed by 
emitting quasars and other AGNs arou nd the 
peak of their activity at z sa 2.5 (e.g., McQii 
|2009 ; Haardt fe Madau 2012 ), later than the 
tion of H T and He T _ It is the traditional vim 
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Vhile plausible, this “two-components” picture for 
coimic reionization is not necessarily correct, and war- 
ra its further investigation for validation. In particular, 
star- forming galaxies at z > 5 can keep the Universe sub¬ 
stantially ionized only if one extrapolates the steep UV 
luminosity function well below the observed limits and 
as mines a globally-averaged absolute escape fraction of 
Lj man-continuum (LyC) radiation into the early IGM, 
fAr. that exceed s 20% (e.g. 


st( in et al. 2012 ; Haardt & Madau 12012 b Such 


al 


da|t 

at 


bad 


ce: 


et 


at 

as$' 


19991; Meiksinl 2005 ; Faucher-Giguere et al. 2009 ). 
oh 


ouwenset al. 2012; Finkel- 


_ eakage 

ues are higher than typically inferred from observa- 
ns of luminous galaxies at z ~ 3 — 4 once contamina- 
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Despite significant efforts and 
examination of hundreds of galaxies, there ex ists only 


] landful of robust detections as of today (e.g., Siana et 
|2015|). 
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flic Planck Collaboration (2015) has reported a new, 
smaller value, r = 0.066 ± 0.016 of the integrated 
re onization optical depth from low-multipoles polar¬ 
ization, lensing, and high-multipole temperature CMB 

a. This corresponds to a sudden reionization event 
+1.7 


z = 


_ 14 and reduces the need for a large LyC 
ikground at very early times. To gether wit h the rc- 
: lt c laim by |Giallongo et al. ( 2015| ) (see also |Glikman 
al.| |201l| ) of a significant population of faint AGNs 
T< z < 6.5, these facts have prompted us to re- 
ess a scenario in which quasars and active galaxies 
y actually dominate the cosmic reionization process 
all epochs, with normal star-forming galaxies making 
a negligible contribution due to their small leakages. 
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Reionization by AGNs 


We explore this intriguing possibility below, assuming a 
(DmAa,^) = (0.3,0.7,0.045) flat cosmology through¬ 
out with H 0 = 70 kms -1 Mpc -1 . 


2. QSO COMOVING EMISSIVITY 

Figure [lj shows the inferrred quasar/AGN comoving 
emissivity at 1 ryd as a function of redshift. Our mod¬ 
eling is based on a limited number of contemporary, 
optically-sel ected AGN samples (see also Khaire & Sri- 
|anand 2015 for a similar compilation). All the surveys 
cited below provide best-fit luminosity function (LF) pa¬ 
rameters, which are then used to integrate the LF down 
to the same relative limiting luminosity, L m i n /L* = 0 . 01 . 
Most of these LFs have faint-end slopes > —1.7, which 
makes the corresponding volume emissivities rather in¬ 
sensitive to the valu e of the adopted limiting luminosity. 
[Schulze et al. ( 2009| ) combined the Sloan Digital Sky Sur¬ 
vey (SDSS) and the Hamburg/ESO survey results into 
a single z = 0 AGN LF covering 4 orders of magnitude 
in luminosity. In the redshift range 0.68 < z < 3 .0. the 


g-band LF of Palanque-Delabrouille et al. (2013) com¬ 
bines SDSS-III and Multiple Mir ror Telescope quasa r 



data with the 2SLAQ sample of Croom_^^_alJ 
The 1 < z < 4 AGN LF by Bongiorno et al. 
again merges SDSS data at the bright end with a faint 
AGN sample from the VIMOS-VLT Deep Survey. The 
high-redshift quasar LF in the Cosmic Evolution Survey 
(COSMOS) in the bi ns 3.1 < z < 3 . 5 and 3.5 < z < 5 has 
been investigated by Masters et al. ( 2012 ), who find a de¬ 
crease in the space density of faint quasars by roughly a 
factor of four from redshift 3 to 4. A significant ly higher 
numb e r of faint AGNs at z ~ 4 is found by Glikman 
|st al. ( 2011 ) in the NO AO Deep Wide-Field Surv e y and 
the Deep Lens Survey, and by Giallongo et al.| ( 2015| ) 
at * = 4 — 6 in the CANDELS GOODS-South field. A 


novel detection criterion is adopted in Giallongo et al 
( [2015 ), whereby high-redshift galaxies are first selected 


in the NIR H band using photometric redshifts, and be¬ 
come AGN candidates if detected in X-rays by Chandra. 
AGN candidates are found to have X-ray luminosities 
and rest-frame UV/X-ray luminosity ratios that are typ¬ 
ical of Seyfert-like and brighter active nuclei. If correct, 
these claims suggest that AGNs may be a more signif¬ 
icant contributor to the ionizing background radiation 
than previously estimated. 

We have converted the integrated optical emis¬ 
sivity inferred from these studies, e\ (in units of 
ergs -1 Mpc -3 Hz -1 ), into a 1 ryd emissivity, £ 912 , us¬ 
ing a power-law spectral energy distribution, eqy? — 


ex( A/912 ) -Quv / esc , with a uv = 0.61 following [Lusso et 
pi] ( |2015| ). We assume an escape fraction of hydrogen- 
ionizing radiation / esc = 1. To assess whether a faint 
AGN population can dominate the cosmic reionization 
process under reasonable physical assumptions, we adopt 
in the following an AGN comoving emissivity of the form 


log £ 912 ( 2 ) = 25.15e 


— 0.0026z 


- 1.5e 


-1.3z 


(1) 


for z < zqso, and zero otherwise. Despite the significant 
scatter in the data points, this function fits reasonably 



ity of luminous quasars inferred by Hopkins et al (2007) 
(see Fig. jl[). It is also higher compared to previous es¬ 
timates at low redshift, a fact that could con tribute to 
solve the ‘/photon underproduction crisis” of Kollmcier 


et al. (2014) (see also Khaire & Srianand 2015). 


9.5 


age of the universe (Gyr) 
2.5 1.5 1.1 


0.8 



-Etti 



square), and 

curve shows t" 


c AGN comoving ionizing emissivity inferred from, 
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asters et al.L( 2012 ) (red pentagons), Cflikman et al, 


(dotted line). See text tor details. 


lallongo et al. (2015) (green - square] 
given in Equation (111 

AGN emissivity of Hopkins et al. (2007]) is shown for comparison 


(2011) (blue 
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The LyC 


3. REIONIZATION HISTORY 

Reionization is achieved when ionizing sources have ra¬ 
diated at least one LyC photon per atom, and the rate of 
LyC photon production is sufficient to balance radiative 
recombinations. Specifically, the time-dependent ioniza¬ 
tion state of the IGM can be modeled sem i-analytically 
by in t egrating the “reionizatio n equations” ( Madau et al. 


1999; Shapiro Sz Giroux 1987) 


^Qhii _ hion,H _ Qmi 

dt (Tin) £rec,H 

^^Helll h'iorijHe ^Helll 


dt 


(™He) 


t 


rec,He 


( 2 ) 

(3) 


for the volume fractions Q of ionized hydrogen and 
doubly-ionized helium. Here, the angle brackets denote 
a volume average, gas densities are expressed in comov¬ 
ing units, t rec is a characteristic recombination timescale, 
and fiion = f dv(e v /hv) is the injection rate density of 
ionizing radiation, i.e. photons between 1 and 4 ryd in 
the case of H i (hi on ,H) and above 4 ryd for He ii (hi on ,He)- 
We do not explicitly follow the transition from neutral 
to singly-ionized helium, as this occurs nearly simultane¬ 
ously to and cannot be readily decoupled from the reion¬ 
ization of hydrogen. 
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The ODEs above assume that the mean free path of 
UV radiation is always much smaller than the horizon 
and that the absorption of photons above 4 ryd is dom¬ 
inated by He n. Because only a negligible amount of 
recombinations occurs in mostly neutral gas, these equa¬ 
tions do not explicitly account for the presence of op¬ 
tically thick absorbers that reduce the mean free path 
of LyC radiation a nd may further delay reionization 
(Bolton et al. [2009 ). They allow, mathematically, for 
values of Q that are > 1, which is physically impossible. 

Following th e results of cosmo logical h ydrodvnamical 
simul ations by Shull et al. ( 2012|) (see also Finlator et al 


p0l2| ), we define the characteristic hydrogen recombina¬ 
tion timescale as 


i ndependent upper li mit of 1 — Qhii < 0.11 at z = 5.9 
( McGreer et al. 2015| ). 

The redshift-dependent fraction of color-selected 
galaxies revealing Lya emission has become a valuable 
constraint on the evolving neutrality of the early IGM. 
We notice here that, in our late reionization model, the 
hydrogen neutral fraction evolves quite rapidly from red- 
shift 6 (1 — Qhii = 0.15) to redshift 8 (1 — Qhii = 0.65). 
This will cause a fast decline of the mean Lya line trans¬ 
missivity of the IGM and may explain the swift drop 
observed in the space density of Lva emitting galaxies 


at 0 > 6 (e.g., Choudhury et al. 2015; Schenker et al. 
2014|). 


trec.H = [(1 + X)( n H)(l + z) 3 (<Xb(Y)) Crr]~ 


(4) 


where as is the case-B radiative recombination rate coef¬ 
ficient, x — Y/[4(l — Y)] = 0.083 allows for the presence 
of photoelectrons from He n (here Y is the primordial 
helium mass fraction), and Crr = 2.9[(1 + z)/ 6] -11 is 
the clumping factor of ionized hydrogen that accounts for 
both density and temperature effects on the average re¬ 
combination rate. Similarly, the recombination timescale 
of doubly ionized helium is 

tree.He = [(l + 2 X ) <n H > (1 + z ) 3 Z(a B (T / Z 2 )) Crr], (5) 

where Z = 2 is the ionic charge and we have assumed that 
H ii and He hi have similar clumping factors. We have 
numerically integrated the reionization equations from 


While a late hydrogen reionization epoch is not a 
unique feature of an AGN-dominated scenario, the early 
reionization of He ii is. Traditionally, the peak in lumi¬ 
nous quasar activity observed at z « 2.5 has been ex¬ 
pected to coincide with the end of helium reionization. 
And while observations of patchy absorption in the He n 
Lya forest at these epochs have been interpreted as sig- 
naling the tail en d of He n reionization at z < 2.7 (e.g., 
Shull et al. [2010 , and references therein), the large fluc¬ 
tuations m the background radiation field above 4 ryd 
predicted from the rarity of bright quasars and the rel¬ 
atively short attenuation lengths of EUV photons make 
the interpretation of the He n data still controversial. In¬ 
deed, rather than complete Gunn-Peterson absorption at 
higher redshifts, recent observation s have revelead high - 
transmission regions out to z = 3.5 ( Worseck et al. 2015| ), 
in conflict with numerical models of He ii reionization 


-QSO 

= 12 o 

nwards, extrapolating the AGN emissivity in driven by luminous quasars ( 

McQuinn et al. 200E; Com- 

Equi 

tion (1 

) to zqso, and assuming an EU V power-law postella et al. 

2013). According to Worseck et al. ( 

20ld), 


and a gas temperature of T = 20,000 K. The integrated 
electron scattering optical depth can be calculated as 

, . , , f z (1 + z') 2 dz' 

t(z) = CCTT(n H ) J - H(~') -^ HI1 ^ + *) + X^Helll], 

( 6 ) 

where c is the speed of light, <jt the Thomson cross sec¬ 
tion, H(z) is the Hubble parameter, and we assumed 

QHell = QHII — QHelll- 

Figure [2| shows the resulting ionization history, quan¬ 
tified by Qhii, QHein, and t ( z ). The shading shows 
the effects of c hanging clumping fac tor [Crr = 9.25 — 
7.211og(l + z), Finlator et al. 2012|| , IGM temperature 
(T = 15, 000 K) and EUV spectral slo pe (a euv = 1.57, 
Telfer et al. 2002 ; Stevans et al. 2014 ). With our de- 
fauit parameters, hydrogen reionization is completed by 
z cz 5.7, helium is doubly ionized by z cz 4.2, and the 
Thomson scattering optical depth is t cz 0.056. The last 
is consistent with the value reported by the Planck col¬ 
laboration at the lcr level. Th e first agrees with an up¬ 
dated analysis by Becker et al. ( 2015 ) of the line-of-sight 
variance in the intergaiactic Lya opacity at 4 < z < 6, 
showing that the data near z = 5.6 — 5.8 require fluc¬ 
tuations in the volume-weighted hydrogen neutral frac¬ 
tion that are higher than expected from density varia¬ 
tions alone. These fluctuations are most likely driven by 
large-scale variations in the mean free path, a signature 
of the trailing edge of the cosmic reionization epoch. As 
shown in Figure |2|, our reionization history is also con¬ 
sistent with the dark pixel fraction observed in the Lya 
and Ly/3 forest of z > 6 quasars, which provides a model- 


sorption demands that the bulk of intergaiactic helium 
was already doubly ionized at z > 3.5 by a population 
of early EUV sources. The high AGN comoving emis- 
sivity present at z > 4 in our model may accomplish 
just that. We find that singly- and doubly-ionized he¬ 
lium contribute about 13% to r, and the He hi fraction 
is already 50% when hydrogen becomes fully reionized at 
redshift 5.7. 

4. DISCUSSION 

Only the most luminous distant quasars can be de¬ 
tected in surveys such as the SDSS, leaving the contri¬ 
bution of faint AGNs to the early ionizing background 
highly uncertain. Recent multiwavelength deep surveys 
have suggested the prese nce of a hitherto unknown pop- 


al 


ulation of faint AGNs (Fiore et al. 2012; Giallongo et 
201E|) and have prompted a reexamination of the role 


played by AGNs in the reionization of the IGM. We have 
expanded on previous studies and assessed a model in 
which the UV radiation responsible for the reionization of 
hydrogen and helium in the universe arises entirely from 
quasars and active galaxies. We have assumed here that 
normal star-forming galaxies will make only a negligible 
contribution to the AGN LyC emissivity given in Equa¬ 
tion ( 0 ), which is true provided their globally-averaged 
escape fraction does not e xceed a few percent or so (e.g., 
Madau fe Dickinson 2014 ). 

Compared to the standard picture widely discussed in 
the literature, the AGN-dominated scenario examined in 
this work completes hydrogen reionization late (z < 6), 
double helium reionization early (z >4), and produces 
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Fig. 2. — Reionization history for our AGN-dominated scenario. Left panel: evolving H II (blue) and He in (magenta) ionized volume 
fractions for our AGN-dominated scenario. Hydrogen in the IGM is fully reionized when Qhii — 1, while helium is doubly reionized when 
OHeiii — 1* The solid lines correspond to our default model parameters, while the shading shows the effects of changing clumping factor, 
IGM temperature, and EUV spectral slope (see text for details). Note the “non-standard” early reionization of helium, Qnein = 1 at 
z > 4. The data point at z = 7 and the upper limit at z — 8 show the constraint on the neutr al hydrogen fraction o f the IGM inferred from 


the redshift-dependent prevalence of Lya emission in the UV spectra of z = 6 — 8 galaxies (Schenker et al. 
z = 5.9 shows the bound on the neutral hydrogen fraction of the IGM inferred from the dark pixel statistics u 


, limit at 

[VlcOreer et alj|2015[ ). Right 

panel: Thomson optical depth to electron scattering, t, integrated over redshift from the present day (green curve). The Planck constraint 
r = 0.066 d= 0.016 is shown as the gray area. 
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Fig. 3.— The hydrogen photoionization rate, Thi, from z = 2 
to z = 5. The green solid line corresponds to our default AGN- 
dominated model, while the shading shows the effects of changing 
the IGM mean free path to ionizing radiation by ±25 percent. 
Red s quares- empirical measurements from the I.vft forest effective 
opacity by Becker & Boltyaa J, 2SXL2)- ~ 
sar- 


with the Planck value at the ler level. It may provide 
an explanation to some otherwise puzzling recent find¬ 
ings, from the rapid decline of the space density of Lya 
emitting galaxies observed at z > 6, to an IGM whose 
temperature is found to increase from redshift 5 to 2 
and where He n appears to be predominantly ionized at 
2 ~ 3.5. It is, of course, a model that is also plagued by a 
number of uncertanties regarding the properties of faint 
AGNs (their space densities and EUV spectra at high- 2 ), 
and of the early IGM (its temperature an d dum piness). 
As already pointed out by Giallongo et al. ( 2015| ), a crit¬ 
ical assumption is the high escape fraction of ionizing 
radiation needed for the global AGN population to dom¬ 
inate reionization. While no discernible continuum edge 
at 912 A is seen in a composite FUV spectrum of 159 
0 < z < 1.5 quasars and active galaxies obtained with 
the C osmic Origins Spectrograph (COS) ( [Stevans et al. 


2014), it is under whether escape fractions of order unity 
are also typical of fainter, higher redshift AGNs. The 
model is also somewhat sensitive to zqscp the “forma¬ 
tion epoch” of the earliest AGNs. Here, we have assumed 
~QSO = 12 and a comoving AGN emissivity that is only a 
factor of 3 smaller at zqso than at redshift 4. In a model 
with zqso = 9, for example, hydrogen reionization would 
only be completed by 2 ~ 5.4. 

As a corroboration of the slowly evolving LyC emissiv¬ 
ity adopted in this work, we have run a mo dified ver¬ 


sion of our radiative transfer code CUBA (Haardt fe 


“r w ^ - -1 - g -~i rr—■— 

quas ar proximity effect ^ Calvcrley et alj | 2011 


SamC usmg tlic _ Madau 2012 ) to estimate the intensity and spectrum of 


the filtered ionizing background - specifically the hydro¬ 
gen photoionization rate Thi - predicted by our AGN- 


a low electron scattering optical depth that is consistent 
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dominated scenario, 


r ffl = 


r, 4^, 

/ av— - a m [v). 

L nl hu 


Here, is the angle-averaged monochromatic intensity, 
h is the Planck constant, and <thi and r'm are the hydro¬ 
gen photoionization cross section and ionizing threshold 
freq uency. We have' computed the background ionizing 

inter sity J using the emissivity in Equation (p]h a EUV 

spectrum with a eU v = 1.7, and a new a piecewise pa¬ 
rameterization of the distribution in redshift and column 
density of intergalactic absorbers that fits the measure- 
ments of the mean free path of 1 ryd photons by Worseck 
et al J (2014). The results, shown in Figure p|, are in for- 
mal agreement with empirical determinations of Thi in 
the in terval 2 < z < 5 based on the Lya opacity of the 
IGM ( Becker fe Bolton 2013| ). Note that, in the tra¬ 
ditional view, the slowly evolving emissivity needed to 
reproduce the Lya opacity data can only be achieved 
by carefully tuning the escape fraction from star-forming 
galaxies to increase rapidly with lookback time, so as to 
compensate for the decline in t he star formation activ¬ 
ity towards early epochs fe.g., lH aardt & Madan| 2012 


poo 

Hi = {rii/n h) / J v <Ji(hv - hvi)dv/{hv), (8) 
for species i = Hi, He 'q or He ii. Here, <7, and Vi are the 
(7) photoionization cross section and threshold frequency for 
the respective species, and n, is their abundance (com¬ 
puted integrating the non-equilibrium rate equations for 
gas at the mean density). The formula above should 
provide the correct mean heating rate of intergalactic 
gas once the background intensity J v is pr operly repro- 


Kuhlen & Faucher-Giguert 2012). 


age of the universe (Gyr) 
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red sh iff 


cessed w hile propagating through the IGM (Puchwein et 
al. 2015|). He ii photoheating exceeds the hydrogen term 


by nearly a factor of three. Measurements of the IGM 
temperature evolution from redshift 5 to 2 derived from 
the Lya forest have been known to be inconsistent with 
the mononotic decrease with redshift expected after th e 


Fig. 4.— Photoheating rates for H I (blue curve). He I (green 
curve), and He II (red curve) in our AGN-dominated model, from 
the present epoch to z = 5, i.e. after the reionization of hydrogen. 
All photoheating rates are expressed per hydrogen atom in units 

of eV |s, 1 , and aiu cumpulud fui gas at l lie iikmii l usniii density— 


completion of hydrogen reionization (Becker et al. 
This discrepancy may be solved by the additiona. 


201l| ). 

heat¬ 


ing provided by an earlier and more extended period of 
He ii reionization. 

The population of faint, high redshift AGNs invoked 
here should lea ve an imprint on the cosmic X-ray back¬ 
ground (XRB). Moretti et al. (2012) derived a value of 
J 2 keV — 2.91^3 x 10 -27 erg cm -2 s -1 Hz -1 sr -1 to the 
2keV unresolved XRB. The expected contribution at 2 
keV from AGNs above redshift z x can be estimated as 


(Haardt & Salvaterra 2015) 


J 2 keV = 


-f 

4?r Jzt 


dz 


(1 +z)H(z) 


£2kev(^)(l + Z)~ 


(9) 


where the specific comoving emissivity at 2 keV is related 
to the LyC emissivity by 

, , D ( 2500 N \“ UV ( 2keV 

£2 keV (2 ) = £ 912 G) X Rjj I - -- 

2 ke V v 9121, V 912 7 V2500A 

( 10 ) 

Here, a ox is the optical-to-X-ray spectral index needed 
to join the emissivity at 2500 A with that at 2 keV, 
and Ru is a correction factor that accounts for the pos¬ 
sible contribution of UV obscured (“Type II”) AGNs 
at 2 > z x to the obse rved XRB. Using a uv = 0.61 


as before, a ox = 1.37 (Lusso et al. 


z x = 5, and Ii\\ = 2 flMerloni et al 

J 2k eV = L64xlO- 27 


erg cm 1 s 


20101) a x = 0.9, 
120141) we obtain 


1 Hz 7 sr 


i.e. a con¬ 
tribution of nearly 60% to the unresolved XRB. We find 
that z > 5 active galaxies can reionize the universe with¬ 
out overproducing the unresolved XRB provided their 
properties (i.e. fraction of obscured objects, optical-to- 
X-ray spectral indices) ar e similar to those of their lower 
redshift counterp arts (cf. |Haardt fe Salvaterra 2015; Di- 


jkstra et al. 2004). 

Finally, we conclude by pointing out that, in order to 
promote further testing of this model against new and 
old observations, we will make the results of our ra¬ 
diative transfer calculations of an AGN-dominated UV 


background freely available for p ublic use at http: / / 


www.ucolick.org/~pmadau/CUBA. 


In our AGN-dominated scenario, the heat input from 
helium reionization will start dominating the thermal 
balance of the IGM earlier than in the standard pic¬ 
ture. In Figure 0 we plot the photoheating rates after 
the reionization of hydrogen for gas at the mean density, 
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